Manipulating atoms in an optical lattice: Fractional fermion number 
and its optical quantum measurement 
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We provide a detailed analysis of our previously proposed scheme [Phys. Rev. Lett. 88, 180401, 
(2002)] to engineer the profile of the hopping amplitudes for atomic gases in a ID optical lattice 
so that the particle number becomes fractional. We consider a constructed system of a dilute two- 
species gas of fermionic atoms where the two components are coupled via a coherent electromagnetic 
field with a topologically nontrivial phase profile. We show both analytically and numerically 
how the resulting atomic Hamiltonian in a prepared dimerized optical lattice with a defect in the 
pattern of alternating hopping amplitudes exhibits a fractional fermion number. In particular, in 
the low-energy limit we demonstrate the equivalence of the atomic Hamiltonian to a relativistic 
Dirac Hamiltonian describing fractionalization in quantum field theory. Expanding on our earlier 
argument [Phys. Rev. Lett. 91, 150404 (2003)] we show how the fractional eigenvalues of the particle 
number operator can be detected via light scattering. In particular, we show how scattering of far-off 
resonant light can convey information about the counting statistics of the atoms in an optical lattice, 
including state-selective atom density profiles and atom number fluctuations. Optical detection 
could provide a truly quantum mechanical measurement of the particle number fractionalization in 
a dilute atomic gas. 

PACS numbers: 03.75.Ss,03.75.Lm,05.30.Pr,11.27.-|-d 



I. INTRODUCTION 



Loading alkali-metal atoms into optical lattices has 
opened up fascinating possibilities to study many-particle 
systems with a wide range of technology to manipulate 
the parameters of the system Hamiltonian as well as the 
quantum states of the atoms. In particular, the flexibil- 
ity of the experimental preparation of ultra-cold atoms 
could be an advantage in realizing atomic systems with 
close analogies to the physical systems, known in other 
disciplines of modern physics, in which the experimental 
evidence or theoretical understanding is limited. 

We previously proposed a scheme to engineer the spa- 
tial profiles of the hopping amplitudes for atoms in an 
optical lattice [l|. This is achieved by superimposing 
coherent electromagnetic (em) field amplitudes that in- 
duce a coupling between different spin states of the atom 
and drive transitions between neighboring lattice sites. 
We considered a particular constructed system of a two- 
species gas of fermionic atoms in a ID lattice, where 
the two components are coupled via an em field with 
a topologically nontrivial phase profile (with topological 
properties similar to a soliton, or a phase kink), so that 
the particle number becomes fractional. In particular, 
we showed that, in the low energy limit, the resulting 
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Hamiltonian of the atomic system has a one-to-one corre- 
spondence to a relativistic Dirac Hamiltonian describing 
fractional fermions in quantum field theory 0, H, I3| ■ This 
is related to fractionalization in the polyacetylenc poly- 
mer systems, where the linearized lattice vibrations are 
coupled to the electrons so that the dynamics of the elec- 
trons is described by a relativistic Dirac equation. We 
also showed that the fractional part of the particle num- 
ber in the atomic system can be accurately controlled by 
modifying the effective detuning of the cm field. 

Moreover, we proposed an optical measurement 
scheme for fractional particle number in an atomic gas 
that could detect directly the fractional eigenvalue and 
determine atom number fluctuations In particu- 

lar, the fractional fermion could be imaged in a phase- 
contrast set-up, while the atom number fluctuations are 
reflected in the fluctuations of the number of photons in 
the light scattered off the lattice. 

In this paper we present a detailed description of how 
fractionalization manifests itself and how it could be op- 
tically detected in the atomic regime, using an optically 
trapped Fermi-Dirac (FD) atomic gas. We show how 
to engineer spatially inhomogeneous hopping amplitude 
profile for the atoms in an optical lattice. We extend our 
analysis in Ref. [l[ of different optical and rf/microwave 
fields to prepare first a dimerized lattice and later to in- 
clude a discontinuity in the hopping amplitude pattern. 
We demonstrate both analytically and numerically how 
the resulting atomic Hamiltonian exhibits a fractional 
fermion number. We consider different observables, ac- 
cessible via light scattering experiments, that could di- 
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rcctly probe the fractional fcrmion number and its fluc- 
tuations. 

Adjusting the hopping amphtudes in our scheme in 
Ref. [l| provided the means of preparing a topologically 
non-trivial atomic ground state. More recently, similar 
techniques to engineer the spatial profile of the hopping 
amplitude between the atoms in adjacent optical lattice 
sites using em transitions were proposed to create effec- 
tive magnetic fields in a 2D lattice @, 01 j as well as non- 
Abelian gauge fields for the atoms in an optical lattice 
[1] . In Ref. [a] a spatially alternating pattern of hopping 
amplitudes along one spatial direction was used to induce 
a non- vanishing phase of particles moving along a closed 
path on the lattice. Since this phase is proportional to the 
enclosed area, the hopping amplitude pattern simulates 
the effects of a magnetic field on charged particles. By 
introducing a level degeneracy in such a set-up allowed to 
consider non-Abelian phases for the atoms moving along 
closed paths [1]. 

An interesting property of our constructed optical lat- 
tice Hamiltonian is that the ultra-cold fermionic atoms 
exhibit long-wavelength dynamics that is equivalent to 
relativistic Dirac fermions with a controllable mass and 
with coupling to a bosonic field. Here the bosonic field is 
represented by the external em field amplitude inducing 
the hopping of the atoms between adjacent lattice sites 
and the mass of the Dirac fermions is proportional to 
the effective detuning of the coupling field from the res- 
onance. The realizations of emerging relativistic Dirac 
quasiparticles in condensed matter systems has recently 
attracted considerable interest, e.g., in graphene sys- 
tems where the 2D hexagonal lattice structure has low- 
energy excitations analogous to those of Dirac fermions 
d, [l0[ [HI, . The atomic optical lattice systems can be 
engineered and controlled to a much higher degree than 
graphene or polymer systems. For instance, the atomic 
gas could be used to demonstrate the effect of the Dirac 
fcrmion mass on the fractional part of the eigenvalue [Tsj 
and the finite-temperature contributions to the fractional 
fermion number in relativistic field theory [3, EB] • The 
optical lattice could provide a laboratory to simulate the 
properties of relativistic Dirac fermions also more gener- 
ally 

Particle number fractionalization is a remarkable phe- 
nomenon in both relativistic quantum field theory and 
condensed matter systems [1, [l^ . Jackiw and Rebbi 0, Q 
showed that for a fermionic field coupled to a bosonic 
field with a topologically nontrivial soliton profile, the 
fermion number can be fractional. The noninteger par- 
ticle number eigenvalues may be understood in terms of 
the deformations of the Dirac sea (or the hole sea) due 
to its interaction with the topologically nontrivial envi- 
ronment. Fractional fermion number has been discussed 
previously in the condensed matter regime in ID conju- 
gated polymers (the SSH model) [1, E The exis- 
tence of fractionally charged excitations in the polymers 
is typically demonstrated indirectly by detecting the re- 
versed spin-charge relation [isj . The fractional quan- 



tum Hall effect (FQHE) can also be explained by invok- 
ing quasiparticles, each with a fraction of an electron's 
charge but the fractionalization mechanism is very 
different from that in the polymers and in the atomic 
gas in this paper. The fluctuations of the tunneling cur- 
rent in low-temperature FQHE regime have been mea- 
sured (20l . [2l| . Interpreting the current shot noise accord- 
ing to the Johnson-Nyquist formula duly suggests that 
the current is carried by the fractional Laughlin quasi- 
particles. Analogous experiments have determined the 
fractional expectation value of the charge in FQHE in 
the Coulomb blockade regime [l^. Recently fractional 
statistics was observed by realiz ing a quasiparticle inter- 
ferometer in the FQHE regime [23j . 

Light scattering can be an efficient probe of quantum 
statistical properties of ultra-cold atoms. It can provide 
information about spatial correlation functions between 
the atoms [l^l and spectral properties about excitations 
[2^. The scattering process could also generate an en- 
tanglement between the photons and the many-particle 
atomic state [2^ In Ref. Q we adapted the tech- 
niques of far off-resonant light scattering in optically thin 
samples [2^ to the atomic gases trapped in optical lat- 
tices. In the measurements of the expectation value of 
the fractional particle number, we make use of the fact 
that the adjacent lattice sites are occupied by differ- 
ent fermionic species. In a phase-contrast imaging, we 
may adjust the light scattered from one species to add 
to the incident beam and the light scattered from the 
other species to subtract. By measuring the intensity 
of the scattered light, one may detect the atom num- 
ber ffuctuations in order to demonstrate that not only 
the expectation value is fractional but also that the fluc- 
tuations become vanishingly small. This is required to 
show that the fractional fermion number represents an 
eigenstate of a particle number operator. The detec- 
tion of both the expectation value and the fluctuations 
would represent a truly quantum mechanical measure- 
ment of fractionalization that has been a long-standing 
challenge in physics. Our proposed measurement tech- 
niques could also be useful in detecting fractional Laugh- 
lin quasi-particles in rapidly rotating atomic gases [28|. 
Although the emphasis of our optical detection analysis 
was on fractionalization, the basic formalism could also 
be used as a more general optical probing technique of 
atom number fluctuations in optical lattices. The optical 
detection of atom number fluctuations in an optical lat- 
tice using an analogous approach and the enhancement 
of the signal by an optical cavity was more recently ad- 
dressed in Ref. [2^, [lol . 

In addition to providing atomic physics technology for 
a potentially direct detection of the fractional fermion 
number and its fluctuations that have so far escaped ex- 
perimental observation, our dilute atomic gas has a possi- 
ble advantage, compared to condensed matter systems, in 
the sense that the atoms form a very clean system: The 
interatomic interactions are weak and well-understood, 
the lattice may be prepared without imperfections, and 
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different parameters of the system Hamiltonian can be 
controlled. 

In Sec. ini we introduce the basic optical lattice sys- 
tem to produce a fractional fcrmion number. In Sec. Ill Al 
we first explain the lattice Hamiltonian and em-induced 
atom hopping from site to site, then discuss the nonuni- 
form hopping amplitudes between adjacent lattice sites 
in a dimerized lattice. Introducing a defect in the reg- 
ular pattern of hopping amplitudes results in a bound 
state and a fractional fermion number, localized around 
the defect; methods to prepare such a defect are also dis- 
cussed. Pertinent mathematical properties of the frac- 
tional fermion number and its fluctuations in a finite lat- 
tice are derived in Sec. Ill Bl An example hopping am- 
plitude profile that generates a fractional number 1 /3 is 
introduced in Sec. lIIII We show that the long- wavelength 
limit of the atomic lattice Hamiltonian is equivalent to 
the relativistic Dirac Hamiltonian in Sec. IIVI We re- 
view some basic properties of the corresponding Jackiw- 
Rebbi relativistic model that provides a simple descrip- 
tion of the fractionalization process. The optical de- 
tection method is explained in Sec. |Vl We first show 
how the fractional fermion number could be imaged in 
a phase-contrast set-up and then how the fiuctuations of 
the fractional observable are reflected in the fluctuations 
of incoherently scattered light. Finally, a few concluding 
remarks are made in Sec. I VII where we also address the 
effects of the harmonic trapping potential. 



II. OPTICAL LATTICE SYSTEM 
A. Basic physical system 

In our scheme to realize particle number fractionaliza- 
tion in atomic gases we consider neutral FD atoms loaded 
in a periodic ID optical lattice. The optical potential is 
induced by means of the ac Stark effect of off-resonant 
laser beams [sTj . 

Atomic systems in ID may be created by confining the 
atoms tightly along the transverse directions. This is 
typically obtained by means of optical lattices or atom 
chips. For instance, in recent experiments ID optical lat- 
tices were created from 3D optical lattices by increasing 
the optical lattice barrier height along two orthogonal di- 
rections [H, [11] . This resulted in a 2D array of decoupled 
ID tubes of atoms where the atoms in each tube expe- 
rienced a periodic optical standing wave potential along 
the axial direction of the tube. If bosonic atoms are con- 
fined sufficiently tightly along the two transverse direc- 
tions, the system approaches the Tonks-Girardeau regime 
where the impenetrable bosons obey the FD statistics 

We analyze an optical lattice system confining two 
fermionic species using an atomic Hubbard Hamiltonian 
[33 |. In a sufficiently deep lattice each site is assumed to 
support one mode function (the Wannier function) that is 
weakly coupled to two nearest neighbor sites. Wc assume 




FIG. 1; The energy diagram of tlie two-species Fermi-Dirac 
gas in an optical lattice. The atoms occupy two different inter- 
nal levels and experience different periodic optical potentials 
shifted by d (on the left). The coupling between the neigh- 
boring lattice sites, representing different spin components, is 
induced by two-photon Raman transitions, or by a superposi- 
tion of the Raman transition and a one-photon microwave or 
rf transition (dotted line). The transition region is denoted 
by the dark shaded area representing the overlap region of 
the atomic wave functions. The energy difference between 
the two components in the overlap region is denoted by u>. 
The two-photon transition is far-detuned by A from an inter- 
mediate atomic level (on the right), and in this illustration 5 
stands for the two-photon detuning. 



that the hopping of the atoms between adjacent lattice 
sites only occurs as a result of driving by coherent em 
fields. This could be realized, e.g., with a FD gas occu- 
pying two internal levels | f) and | j) that are coherently 
coupled via an cm-induced transition, if the two species 
experience optical potentials which arc shifted relative to 
each other by d, where 2d denotes the lattice period for 
each state | t) and | |); see Fig. [T] The coupling could 
be a far off-resonant optical Raman transition via an in- 
termediate atomic level, a microwave, or a rf transition. 
The optical lattice potential with alternating spin states 
in alternating sites may be realized, e.g., by combining 
(7"*" and a~ polarized lasers [ssj . or when the light beam 
is blue-detuned from an internal transition of the atoms 
in level | t) and red-detuned from a transition in level 
I J,). A simple example ID lattice potential in that case 
is V^{x) = Fosin2(7rx/2d), and Vi{x) = -Vosm'^{nx/2d). 

In our example set-up of Fig. [U the neighboring lat- 
tice sites represent atoms in different internal levels and 
are separated by a distance d. Here d = A/[4 sin(0/2)] is 
obtained by two laser beams, with the wavelength A, in- 
tersecting at an angle 6. The lattice spacing can be easily 
modified by changing the angle between the lasers. We 
may write the annihilation operators for fermionic atoms 
as C2fc+i,| and C2k.i at odd and even numbered lattice 
sites, occupying the internal levels t and |, respectively. 
In the following we suppress the explicit reference to the 
different internal levels and write the annihilation oper- 
ator for the atoms in the fcth site simply as c^.. The 
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Hamiltonian for the atomic system then reads 

7^ = H [^kclck - (KfcC^iCfe + 4c3;Cfe+l)J . (1) 
k 

The em-induced coupHng between the two internal states 
between the /cth and the {k + l)th site is described by 

nk= j d\ra{^-Vk)n{v)ijj_„{Y~Vk+i), (2) 

and Ck stands for the effective detuning and/or exter- 
nal potential variation between different sites. The em- 
coupling terms are the analogues of the hopping terms in 
the corresponding polymer Hamiltonian [3|. The mode 
functions of the individual lattice sites (Wannier func- 
tions) are denoted by '(/'^(r ~ r^), with a =t,i and 
— fT =i,t, respectively. We assume that the em cou- 
pling between the internal levels r2(r) is the only transi- 
tion mechanism for the atoms between neighboring lat- 
tice sites, and therefore ignore the direct tunneling. ri(r) 
could, for instance, be an effective two-photon Rabi fre- 
quency at the position r. In this respect our optical- 
lattice scheme is closely related to the microtrap scheme 
discussed recently as a method of transporting atoms for 
the purposes of quantum information processing [3^ . We 
neglect the s-wave scattering between the two FD species. 
This process can be made weak by having a small spa- 
tial overlap between the mode functions of the adjacent 
lattice sites. 

Even though the optical lattice may be part of a larger 
trap, which could generate interesting physics in its own 
right, we simphfy by setting = [—1^5/2 to be con- 
stant along the lattice separately for the even (|) and for 
the odd (t) sites that are split by the energy difference 
5. We mostly concentrate on the simplest case (5 = 0. 
The fractional fermion number arises from certain types 
of defects in the couplings k^. In this paper we illus- 
trate by concentrating on one such set-up; we consider 
a dimerized lattice generated by the coupling matrix ele- 
ment that alternates from site to site between two values 
a + ^ and a ~ fi, except that at the center of the lattice 
there is a defect such that the same coupling matrix ele- 
ment appears twice. In the following, we use the notation 
introduced in Ref. 0, [s^] • 

1. Dimerized optical lattice 

In this section wc introduce an alternating pattern in 
the atomic hopping amplitudes. This establishes a dimer- 
ized optical lattice. If we also add a defect in the pattern 
of alternating coupling matrix elements, this results in a 
fractional fermion number, localized around the defect. 
This type of schemes are the subject of the following two 
sections. 

We consider the coupling frequency ^{r) to be a phase- 
coherent superposition of two field amplitudes: 

r!(r) = V(r)+Z^(r)sin(^) . (3) 



sin(7rx/d) 




a+fi a-n a+ij a-fi a+jj a-fi 

FIG. 2: Dimerized optical lattice with spatially alternating 
hopping amplitude for atoms between adjacent lattice sites. 
The hopping is induced by electromagnetic field (as in Fig. [l| 
whose amplitude varies according to Eqs. ([Sjl and (|6}. Here 
ak and fik are assumed to be constant along the lattice. The 
coefficient in the front of /x is determined by the value of 
sin(7ra;/d) at the middle between the lattice sites where the 
Wannier functions of the two lattice sites overlap. 

We can assume that the matrix elements 

ak= J d^'rrAr - rfe)V(r)V^_,(r - rfc+i), (4) 
/I, = y d^rAr - Vk)Uiv)^-Ar - rk+i) , (5) 

do not change sign over the length of the lattice. The 
purpose of the standing wave in Eq. ([3|) is to introduce 
a spatially alternating sign for the hopping amplitude 
along the lattice. 

We assume that sin(7rx/c?) in Eq. ([3]) is approximately 
constant over the small spatial overlap area of two neigh- 
boring lattice site atom wavefunctions. One of the lattice 
sites has the potential minimum at a; = and its two 
nearest neighbor sites minima at a; = id. The spatial 
overlap of the atomic wavefunctions centered at the sites 
X = and x = ±d has the maximum at a; = ±d/2. The 
sin(7r2;/c?) term in the hopping amplitude in the overlap 
region at x = d/2 is equal to one and at a; = —d/2 it is 
equal to — 1 . With similar arguments it is easy to see that 
the hopping amplitude alternates between the values ±1 
at each overlap area and we obtain (see Fig. [2]) 

Kfc ~ flfc + fik ■ (6) 

In this approximation, we take and real and in the 
actual calculations assume them to be constant along the 
lattice, Ofc = a and /.j^ = /i. 

The field amplitude, Eq. ([3]), required to generate the 
hopping amplitude, Eq. ([6]), represents a phase-coherent 
superposition of a standing wave and a field amplitude 
with a uniform phase profile. The field superposition 
could be produced, e.g., by using optical holograms. We 
will address such techniques in more detail in the next 
section, when we consider a more complicated field pro- 
file needed to prepare fractional fermions. Due to the 
relatively simple form of Eq. ([3]) , the field configuration 
could be generated by superposing an optical standing 
wave and a rf, microwave, or an optical field V. The two 
fields both need to couple the states | 1) and | j), or to 
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form part of a multi-level, multi-photon transition that 
couples the two states together. 

The standing wave in the hopping amplitude has the 
period 2d when the period of the amplitude of the lattice 
lasers forming the optical potential is equal to Ad. The 
different periodicity may be obtained by using a different 
angle 0' between two counter-propagating lasers, so that 
d = A'/[2 sin(6''/2)], where A' is the wavelength of the 
sin(7ra;/d) hopping field. 

Alternatively, the sin(7ra;/c?) hopping term may be pre- 
pared by using the same wavelength A' = A and inter- 
section angle 0' = as in the optical lattice potential 
with a two-photon optical Raman transition (Fig. [1]). 
The strength of an off-resonant two-photon Rabi fre- 
quency in the limit of large detuning, A, from the in- 
termediate state is f2 cx 7?.i7?.2/A, where TZi denote 
the Rabi frequencies in the individual transitions [25j . 
For two standing-wave one-photon couplings displaced 
from one another by d, the two-photon Rabi frequency is 
Q. (X sin(7ra;/2(i) cos{'kx /2d) (x sin(7ra;/d), so, even though 
the period of the lattice laser amplitude is 4(i, the effec- 
tive field amplitude of the two-photon transition has the 
desired period 2d. 

In order to understand the effect of dimerization, we 
follow and consider the Hamiltonian ([T]) for the uni- 
form lattice Efe = with the coupling © in the case of 
periodic boundary conditions. The adjacent lattice sites 
are separated by d and we can obtain the reduced zone 
representation for ~T:/2d < p < TT/2d hy transforming 
the lattice operators according to 

Cp+ = -^^e^''^-'^(-lfcfc. (8) 

In terms of the valence and the conduction band opera- 
tors of the undimerized system, Cp_ and Cp+, the Hamil- 
tonian ([T]) reads: 

-J- = ''^^[^p{'^p+'^P+ ~ Cp_Cp_) 

p 

+ Ap(c^+Cp_ + 4_Cp+)] , (9) 

with 

Ap = 2yU sin (pd) , ep = 2a cos {pd) . (10) 

This can be diagonalized by the Bogoliubov transforma- 
tion 

bp- = apCp_ - PpCp+, (11) 
bp+ = apCp+ + (3pCp- , (12) 

with the appropriate choice of the coefficients ap and /3p, 
so that 

|=^£;p(6t^&p+-6j_&p_). (13) 
p 



The quasiparticle energies are given by 

Ep = + A2 , (14) 

and the spectrum exhibits an energy gap equal to 4/i|/i| 
at the reduced zone boundary p — ±7r/2d, due to the 
dimerization field. We may define the atomic correlation 
length according to Eq. analogously to the BCS 

theory: 



e'iPF) 




ad 


A{pf) 







(15) 



where pf = /2d is the Fermi momentum. 

2. Defect in the pattern of alternating hopping amplitudes 

In this section we construct a field amplitude that 
breaks the symmetry of the dimerized lattice formed by 
the alternating pattern of the hopping amplitudes. This 
is obtained by generating a defect at the center of the lat- 
tice, so that the same coupling term between the lattice 
sites appears twice in a row. 

For an even number of sites, or in the limit of an infinite 
lattice, the two hopping amplitude profiles 

Kk = a±{-lfn, (16) 

yield ground states with equal energies. We may join 
these two field configurations at the center of the lat- 
tice by synthesizing a defect in the alternation pattern in 
such a way that the left side of the lattice represents the 
lower sign of the hopping term (flBl) and the right side 
of the lattice the upper sign; see Fig. [31 The em field 
required to induce the desired coupling may be written 
in the following form: 

= V(r) + U{v)v{x) sin (^) , (17) 

where the field <p(x), satisfying ip{h) ~ —ip{—h) = 1 for 
all h ^ d, exhibits a spatial profile of a soliton, or a 
phase kink. Here '^{x) has a phase jump of tt at .t = 0, 
representing a topological phase singularity. The corre- 
sponding hopping amplitude reads 

Kk ^ a + ipk{-l)'' ^l , (18) 

where ipk = v{xk)\ see Fig. [3] 

This type of a coupling represents a superposition of a 
field with a constant phase profile V (producing a) and 
the field Lp{x)sva.{'Kx/d). The entire field configuration 
P?)) could be prepared optically by synthesizing multi- 
beam superposition states using diffractive optical com- 
ponents. In particular, computer-generated holograms 
and spatial light modulators can act as optical phase 
holograms to shape laser fields in order to generate a 
sharp phase kink (p{x) ~ a;/|a;|, combined with the sine 
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-sin(;rx/d) 




sin(;rx/d) 




a+ij I a+ij a-fi a+ju 



FIG. 3: Defect in the pattern of alternating hopping ampli- 
tudes. The coefficient in the front of fj. is determined by the 
product ^p{x) sm{nx/d) at the middle between the lattice sites 
where the Wannier functions of the two lattice sites overlap. 
The soliton, or a kink, profile of ip{x) generates a defect in the 
hopping amplitude pattern, so that the same coupling matrix 
element appears twice. The defect is marked by a dashed 
(red) vertical line. Here ip{x) ~ a;/|2;|. 



wave and the field V. The holograms are typically cal- 
culated by specific algorithms that optimize the phase 
and/or intensity profiles of the field according to the 
merit function criteria. Optical holography has been ex- 
perimentally used to prepare phase singularities on light 
beams jssj and has been proposed as an efficient tool 
for phase imprinting top olog ical defects in atomic Bose- 
Einstein condensates [33, 13 ■ 

Alternatively, one may generate the effective coupling 
field, Eq. ([T7|) . by superposing two different transition 
paths between the states | t) and | |). One of the tran- 
sition paths, that may involve a multi-photon transition 
via intermediate levels or a single-photon transition with 
no intermediate levels, produces an effective Rabi fre- 
quency with a uniform phase profile, representing the 
field V in Eq. and a in Eq. ([TH]). The other transi- 
tion path generates the field ip{x) smijrx/d). This could 
be obtained, e.g., with a two-photon transition via an 
intermediate level, so that the first transition is driven 
by (p{x) and this is followed by the second one driven by 
sin(7ra;/d); see Fig.|H 

The effective (p{x) field might be produced either by 
making use of the (rf or microwave) transition between 
the spin states, or by means of an optical transition other 
than the one used to produce the sin(7rx/d) standing 
wave. In our studies of fractionalization, the precise form 
of the field (p{x) is not very crucial, as the relevant physics 
depends on its asymptotic, topological behavior. A phase 
profile with topological properties similar to 1^9(2;) could 
be prepared, e.g., by means of a standing em wave sin((3'a;) 
satisfying q <C ir/d. The period of sin(ga;) needs to be 
chosen long so that along the lattice there is only one 
node point located near the center of the lattice. Such 
a ID standing wave could be obtained using microwaves 
or, alternatively, optical fields if the intersection angle 
between the two laser beams could be chosen such that 
sin(6'/2) ^ X/2d, where A is the laser wavelength. In 
Ref. [l| we proposed also other potential experimental 
schemes to prepare the hopping field based on em field 




FIG. 4: Examples of the field configurations that are not 
based on an optical hologram and that can generate the field 
amplitude in Eq. (|17|l . Two different transition paths syn- 
thesize the superposition between V and fix) sm{nx/d). The 
field profile <p(a;) sin(7ra;/d) is produced by two successive tran- 
sitions, so that the effective Rabi frequency is the product of 
ip{x) and sin(7ra:/d). On top, a far-detuned microwave or rf 
field, generating ip{x), is combined with an optical Raman 
two-photon transition via a far-detuned electronically excited 
intermediate level, where one of the transitions is driven by 
the standing wave sin(7ra;/d). The effective Rabi frequency for 
the corresponding three-photon transition is proportional to 
cx (fi{x)sm{nx/d). On bottom, an optical Raman two-photon 
transition where one transition is driven by ip{x) and the other 
one by sin(7ra;/d). 



couplings. 



3. Alternative scheme based on field gradients 

The scheme introduced in Ref. [l[ to prepare a nonuni- 
form hopping amplitude pattern for the atoms in an op- 
tical lattice has an alternative proposed realization by 
Jaksch and ZoUer Q. In Ref. @ a spatially alternat- 
ing pattern of hopping amplitudes was generated in a 
2D lattice along one spatial direction in order to induce 
a nonvanishing phase on atoms moving along a closed 
path in the lattice, therefore simulating the effect of a 
magnetic field on electrons. This approach uses field gra- 
dients combined with an em-induced coupling. In the 
following we adapt the basic scheme of Ref. Q to our 
fractionalization study. 

As in the previous section, we assume that the hopping 
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of the atoms between adjacent lattice sites is driven by 
em fields. To be specific, let us think of two-photon tran- 
sitions in the Raman configuration with a far-off resonant 
intermediate state. This setup is functionaUy equivalent 
to transitions in an effective two level system with the 
transition frequency lu equal to the difference of the en- 
ergies of the states "f and I, driven by a field whose ef- 
fective frequency is the frequency difference between the 
two laser frequencies. In Fig.[5Ka) we draw the scheme of 
transitions from site to site. The horizontal lines repre- 
sent the energies of the levels t and |, alternating between 
adjacent sites, and resonant transitions between the lev- 
els at the effective frequency of the driving field uj are 
marked with arrows. In this basic lattice the site-to-site 
couplings are the same. 

Next assume that in some manner a potential has been 
added to the atoms that grows linearly along the lattice. 
Correspondingly, the Hamiltonian picks up an additional 
term 

For notational simplicity the lattice sites are numbered 
from —Nh to Nh- Wc specifically pick £^ = hoj, and imag- 
ine tuning the fields driving the transitions to the ef- 
fective frequency 2uj. As demonstrated in Fig. [5]^b), in 
this way every second transition is on resonance. On the 
other hand, every second transition now has the effective 
resonance frequency 0, and the field with the frequency 
2lo assumedly becomes ineffective at transferring atoms 
in such transitions. To compensate, let us add light fields 
to make the effective frequency to drive the transition 
between the sites that the potential Eq. ([TO]) decoupled, 
as in Fig.[S](c). Resonant couplings between all adjacent 
sites have now been established. To prepare for the final 
step of the argument, in Fig.[5lc) we have also indicated 
a few possible transitions without a resonant final level. 
As we denote symbolically by using solid-line and dashed- 
line arrows for the driving fields, the coupling coefficients 
for the alternating transitions need not be equal. At this 
stage we have a dimerized lattice. 

In the last step we modify the added potential to make 
the corresponding force change sign at the center of the 
lattice, so that 

Hs=(. (20) 

replaces Eq. (fTO|) . Figure. [Hd) shows both the level 
scheme for the newly adjusted lattice potential, and the 
possible transitions. This time around two solid- arrow 
transition in a row are on resonance at the center, so that 
we have established a defect in the couplings precisely as 
needed to prepare a half-fermion. 

The technical issue with these ideas is to set up 
the piecewise linear potentials. A potential of the 
type Eq. could be generated, say, with a magnetic 



(a) 




-2-1012 
j 

FIG. 5: A scheme for a dimerized lattice based on the in- 
terplay between internal and center-of-mass energies of the 
atoms, (a) The energies of the relevant internal states (t and 
I) of an atom at lattice sites j — —2, ... ,2 are marked by 
horizontal bars, with resonant transitions between the states 
displayed by solid arrows, (b) After a constant energy differ- 
ence between adjacent lattice sites is added to the atoms, for 
a given set of driving fields only half of the transitions can be 
on resonance, (c) Another set of lasers frequencies is added, 
which render also the remaining transitions resonant (dashed 
arrows) . A dimerized lattice has been prepared. Several tran- 
sitions without a resonant final state are are also indicated, 
(d) In a variation of the scheme in part (b), a constant en- 
ergy difference is added between every state pair to the left of 
j — 0, and the negative of the same energy difference between 
state pairs to the right of j — 0. The same coupling as in part 
(c) now realize a dimerized lattice with a defect. 



field ramp that shifts Zccman levels, using ac stark shifts 
due to added off-rcsonant laser or microwave fields, in 
the case of polar molecules static electric fields, or even 
by accelerating the optical lattice. The reversal of the 
force at the center of the lattice in Eq. (PO)) could, e.g., 
in the case of off-resonant lasers be prepared using opti- 
cal holograms or spatial light modulators. The required 
field gradients may also be generated by shining a laser 
through an absorption plate with a linearly changing op- 
tical thickness. 
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4- Remarks 

We have demonstrated how em field-induced liopping 
of atoms between adjacent lattice sites can be a very 
useful tool in engineering a rich variety of atomic lattice 
Hamiltonians. A nonuniform field amplitude, with a pe- 
riodic spatial profile, is sufficient to prepare a dimerized 
lattice. One may also introduce a defect in the regu- 
lar pattern of hopping amplitudes. We will next show 
that the hopping amplitude profile with such a defect in 
the pattern of alternating coupling matrix elements re- 
sults in a bound state and a fractional fermion number, 
localized around the defect. Such a coupling frequency 
also converts the atomic lattice Hamiltonian ^ in the 
continuum limit into the Dirac Hamiltonian exhibiting 
fractional particle number in rclativistic field theory. 



B. Fractional fermion in a finite lattice 

In this Section we consider the Hamiltonian ^ in the 
presence of the hopping field with a defect, Eq. (fTS)) . 
along the lines of Ref. [i^l , for a finite lattice and with- 
out a continuum limit. We use hard- wall boundary con- 
ditions, i.e., the lattice simply terminates at the edges as 
if the first sites past the ones included in the lattice were 
forced to be empty. We demonstrate how a lump with 
the particle number expectation value of a half emerges 
on top of the a homogeneous background fermion den- 
sity. We also discuss numerical techniques used to study 
the fractional fermion, both as it comes to the expecta- 
tion value and the fluctuations of the atom number, and 
present representative examples. 



1. Elementary half-fermion 
Suppose we can find fermion operators 7„ of the form 

(21) 



k k 



such that the Hamiltonian ([T]) becomes 
H 



J2 ^p^l^p ■ 



(22) 



then we have diagonalized the Hamiltonian. Any one of 
the operators 7^, call it 7, should then satisfy 



H 



1, 



= 0^7. 



(23) 



It follows readily from Eqs. pT|) -(|23 p that the column 
vectors flfc of the matrix Ukp are the orthonormal eigen- 
vectors of the eigenvalue problem 



which also yields the corresponding eigenvalues LOp. As 
in Eq. (jM]) . the couplings Kk may be assumed real with- 
out a loss of generality, and so we take the matrix U to 
be real and orthogonal 



(25) 



We consider a sharp defect t p(x ) ~ x/\x\ in the hopping 
amplitude profile. As in Ref. [42l |. we take the number of 
lattice sites to be Ns = 2Nh -I- 1 = 4ri, -I- 1 , where n itself 
is an integer, and number the sites with integers ranging 
from —Nh to Nh- For illustration, pick n = 2, use an x 
to denote a lattice site, and ± the couplings a ± /x, then 
our lattice with the couplings reads 

X — x-fx — x-l-x + x — x + x — X (26) 

The defect in the middle, two consecutive plus signs, is 
the crux of the matter. 

We assume here a uniform lattice with the effective en- 
ergy splitting between the adjacent sites 6 ~ 0, so that 
e^, = 0. It is easy to see from the structure of Eq. ([M]) 
that if uj is an eigenvalue, then so is —uj (the conjugation 
symmetry); and the eigenvectors transform into one an- 
other by inverting the sign of every second component. 
We will label the eigenvectors as —Nh ■ ■ ■ Nh in ascend- 
ing order of frequency, and assign the labels ±p to such 
± pair of states. Correspondingly, the transformation 
matrix U satisfies 



\Ukp\ = \U,^^p\ 



(27) 



But under our assumptions, the number of eigenvalues 
and eigenstates is odd. The ± symmetry implies that an 
odd number of the eigenvalues must equal zero. Except 
for special values of the couplings a and ^, there is one 
zero eigenvalue. We call the corresponding eigenstate the 
zero state. Provided a and fi have the opposite signs and 
\a\ > 1^1, all odd components in the zero state equal zero 
and the even components are of the form (37j 



a + /i 
a — /i 



\k\/2 



(28) 



The normalizable zero state becomes the narrower, the 
closer in absolute value a and n are. For a broadly dis- 
tributed bound mode, \a\ 2> |a<|, we obtain from Eq. ([28|l 
the limit 



^c. (-1)1^-1/2 cxp(- 

Xo 



\xk ~ xo\/0 



(29) 



where ^ is the atomic correlation length calculated in 
Eq. (fT5|) . Since the size of the zero state depends on the 
relative strength of the superposed electromagnetic fields 
via the correlation length, it could be varied experimen- 
tally. This is different from the polymer case, where the 
size of the bound state is fixed. 

Suppose next that the system is at zero temperature, 
and contains N 



(efe - uj)ak - KkUk+i - Kk-ittk-i = , 



'f = Nh + 1 fcrmions. The exact eigen- 
(24) states p at the ground state arc then filled up to zero state 
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and empty at higher energies, with occupation numbers 
np= or 1, where 

np = illlp) ■ (30) 

The number operator for the fermions at site k corre- 
spondingly reads from Eqs. (PT|) and 

= X! ^kpUkqlllq , (31) 
pq 

SO that the expectation value of the fcrmion number at 
the site k is 

-1 

(cK)= E \U'^p\'^ E \Ukp\' + \Ukof 

p=-Nh p=-Nh 
-1 Nh 

= h E M + 1 E iC^^pl' + lt^fcop. (32) 

Here we have explicitly separated the zero state and used 
the symmetry \Ukp\ = |J7i;__p|. The matrix U^p defines 
a complete basis and we can use the corresponding or- 
thogonality of U by combining the first two terms and 
one-half of the third term on the last line of Eq. ([5^ : 

{Ack) = h E \Uup? + \\Uk,? = \ + h\U'^,?. (33) 

By virtue of the same orthogonality, localized with the 
zero state there is a lump with | l^fcoP = \ fermions 
on top of a uniform background of half a fermion per site. 

This lump is the celebrated half of a fermion 0, [1, 
Q. It is the result of adding a fcrmion to the zero state 
which has a nonvanishing occupation probability only for 
every second site. Therefore the half- fermion rising from 
the constant background resides on every second lattice 
site, here on the even-numbered sites. For the remaining 
lattice sites the mean occupation number is exactly a 
half. 

As a result of the defect the normalizable zero state 
is localized around the phase kink. This zero state cre- 
ates a fractional deficit of states in the valence and the 
conduction bands, as can be observed in Eq. ([5^ . In 
the presence of the conjugation symmetry, the density of 
states is a symmetric function of the energy, and both 
bands have locally a deficit of one-half of a state. 

2. Fermion numbers and their fluctuations 

To demonstrate the precise meaning of the fractional 
fermion, we define two smoothed fermion number oper- 
ators using an envelope function oik that covers the zero 
state around the phase kink; 

^ = E "'=^1^'^' ^ = E (4cfe - \) ■ (34) 

fc k 



The idea is that an experiment with a limited resolution 
may be expected to address many lattice sites at once. 
Alhtough the numbers ak are real in all of our worked- 
out examples, for maximum flexibility we allow them to 
be complex. The difference between the operators N and 
N is that in the latter we have subtracted the constant 
background fermion number ^ at each lattice site. 

Let us so far assume that the numbers ak vary lit- 
tle from site to site, and that their values in a region 
that covers the half-fermion lump is well approximated 
by unity. We then have from Eq. ([55)1 

{N) = \Y.ak\Uko? (35) 

k 

which just reiterates the observation that the lump con- 
tains, on the average, half of a fermion above the uniform 
background. 

There is more to a fractional fermion number, how- 
ever, than the average density. So far we have only dealt 
with the expectation values of the atom numbers. How- 
ever, the fluctuations can be small as well, indicating 
that the half-fermion represents an eigenstate of the op- 
erator N [i^, [4^. Since the fluctuations are the same 
for the operators N and iV, we focus on the former. The 
expectation value of N is 

{N) = Y.^k{clck) 

k 

= ^^akUkpUkqilllq) 
k 

= ^akUlpUp, (36) 
fc 

and for the square wc have 

= E^'^^'^^fc'^'^^;'^') 

kl 

= E (^k^lUkpUkqUlrUls{-lllqllls) 
klpqrs 

= ^ alai [Ul^UfqUpTlq + UkpUkqUlpUlq{np - UpUg)] 
klpq 

= \{N)\^ + {ANf, (37) 
where the part characterizing the fluctuations is 
{ANf ^{N^^N) - |(iV)|2 

= ^alaiUkpUkqUipUiqnp{l - Ug) . (38) 

klpq 

Although wc always use the zero-temperature Fermi sea 
with the occupation numbers Up = or 1 for the jp 
fermions in our explicit examples, this assumption has 
not been used in the analysis of the fluctuations. With 
the appropriate thermal occupation numbers Up, expres- 
sions such as Eq. ([55]) apply also at finite temperature 
0,[ll. 
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FIG. 6; Coherent atom number fluctuations AA^ under a 
Gaussian (solid line) and alternating-sign Gaussian (dashed 
line) envelopes, and q^, and the corresponding incoherent 
atom number fluctuations {AN)i (dotted line) as a function 
of the width of the Gaussian w. The lattice parameters are 
Af, = 1025, Nf = 513, and n = -0.1 a. 

For illustration we use the Gaussian shape, 



We take a lattice with Ns = 1025 sites, pick the param- 
eters n = —0.1a, put in Nf — 513 fcrmions so that the 
zero state is the last filled stated, and find the rms fluc- 
tuations of the fermion number AA'^ as a function of the 
width of the weight function w. The result is shown as 
the solid line on a log-log plot in Fig. [S] The notch around 
w = 1 indicates that at this point the weight factors ak 
start to span several lattice sites. Another break in the 
curve is seen at about w = 10 when the weight function 
covers the whole zero state. Thereafter the fluctuations 
behave as AiV cx w~^^'^. The fermion number N under 
the weight function becomes more sharply deflned as the 
region for averaging grows broader. Finally, at w ^ 500, 
the weights ak effectively cover the entire lattice. The 
fluctuations then decrease even faster with increasing w, 
as is appropriate for the flxed fermion number in the lat- 
tice as a whole. 

Let us next suppose that the weights ak can be ar- 
ranged to alternate in sign from site to site. Specifically, 
we use 

at = (-l)'=e-('^'/™)' . (40) 

Since the half-fcrmion is confined entirely to the even- 
numbered sites, this weight has the effect that in the limit 
w oo the constant background is automatically sub- 
tracted as far as mean atom numbers go; (N) = (N) ~ i. 
Regarding fluctuations, we again use Eq. p8|l . The re- 
sulting AiV is plotted in Fig. [S] as a function of the width 
w as the dashed line for the same problem parameters 
that were used for the single-sign Gaussian a^. Once 



the envelope covers several sites, the fluctuations for the 
models a^ and go separate ways, and in the latter 
case actually increase with increasing width w until the 
whole lattice is covered. 

Our flnal case of fluctuations is the case of simply 
adding the particle number fluctuation at each site with 
the weights jafcp, as if the fluctuations were independent. 
We thus have the "incoherent" fluctuations, 

(AiV)? =^ K|2 ((4c0 - (cK)(cU)) 
fc 

= 5]|afcpt/^p[/2^np(l-n,). (41) 

kpq 

This is what becomes of the "coherent" fluctuations 
of Eq. ([55]) if only the diagonal terms with k = I are kept. 
A physical motivation for such a rule could be that the 
coefficients ak are uncorrclatcd random numbers (with 
zero mean). 

The incoherent fluctuations (AA^)^ are plotted in Fig. [5] 
as a dotted line along with the coherent fluctuations for 
the two envelopes and q;J|^; which envelope is used in 
the incoherent case does not matter. The magnitude of 
the incoherent fluctuations is between the magnitudes of 
the fluctuations for the two coherent cases. As soon as 
the widths of the Gaussians are sufficient to cover several 
lattice sites, the coherent fluctuations for the alternating- 
sing Gaussian are about a factor of times the 
incoherent fluctuations. 

In the standard half-integer fermion number argument 
one uses a smooth envelope such as a^ and subtracts a 
non-fluctuating neutralizing background of ^ charge per 
lattice site, whereupon (TV) i and AA^ with an 
increasing envelope width w. The intermediate regime 
that occurs once the zero state is covered is the crux of the 
matter. Not only does the expectation value of fermion 
number TV equal i, but the fluctuations are small and TV 
therefore has an eigenvalue ^ . 

A perusal of the dotted line for the incoherent fluctu- 
ations in Fig. [6] qualitatively suggests that the fluctua- 
tions of the fermion number at each individual site are 
comparable to what one would expect if each site sim- 
ply had the fermion number of either or 1 with the 
equal probabilities of i. On the other hand, the atom 
number fluctuations appear to be anticorrelated between 
adjacent sites. Thus, if one sums up the atom numbers 
with a weight that has a sign alternating from site to site 
(dashed line), the resulting fluctuations in the summed 
atom number are larger than they would be for uncor- 
rclatcd fluctuations. Also, if one sums over the atom 
numbers with a weight that varies little from site to site 
(solid line) , such anticorrelated fluctuations tend to can- 
cel in the sum. In short, anticorrelations of atom num- 
ber fluctuations between neighboring sites undoubtedly 
contribute to the impression of a sharp eigenvalue for a 
smoothly weighted sum of occupation number operators. 
The role of longer-range correlations (next-nearest sites, 
etc.) deserves further study, but it will not be undertaken 
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here. 



3. Variations of the half-fermion 

The main point to emphasize about the derivation 
of Eq. ([33|) is that the fractional particle number is the 
property of the system as a whole. Adding the last 
fermion that occupies the zero state adds one particle, 
but makes a soliton containing half of a fermion. This 
means that before the addition there was a hole in the 
density distribution with half of a fermion missing, an 
observation easily verified by a calculation similar to the 
one that lead to Eq. ([55]) . The hole clearly depends on 
the shapes of the wave functions and on the occupation 
numbers of the non-zero states, hence so does the half- 
integer fermion number. 

Fractionalization is a robust phenomenon. Something 
akin to a localized zero state occurs as soon as the reg- 
ular alternation of the couplings between adjacent states 
gets out of rhythm around a defect. In particular, the 
defect does not have to be confined to one lattice site; 
this should make the experiments easier. Moreover, a 
spike with about half a fermion standing above the back- 
ground may be seen when the particle number deviates 
from {Ns -I- 1)/2 by, say, ten per cent; and a hole of about 
half a fermion may similarly persist even if the particle 
number is well under (Ns + l)/2. 

The total number of atoms, of course, must remain 
an integer, and the optical lattice has a finite size. In 
our examples with iV^ = 4n -I- 1 sites and Nf = 2n + 1 
fermions, the background occupation of ^ makes 2n + ^ 
fermions, and the half of the fermion in the peak adds 
up to make the correct total fermion number 2n + I. If 
the number of sites is even, the system compensates by 
putting another soliton at the edge of the lattice, simi- 
lar to Shockley edge states If the number of sites 
is even and the number of fermions precisely half of it, 
numerically, all solitons may vanish. The zero state is 
then doubly degenerate, and, as we have not tried to as- 
sert any control over the state in the degenerate subspace 
that gets the last fermion, the numerics decides for us. 

The restriction \a\ > and a and fi exhibiting the 
opposite signs could be relaxed to read |a -t- ^| < \a — 
fi\. If this condition is violated, the primary soliton, the 
counterpart of the half-fermion in our examples, emerges 
at the edge of the lattice, and the possible extra soliton 
at the center. 



III. FERMION NUMBER 1/3 

We may also generalize the half-fermion concept to 
fractionalization to values 1/n, for any integer n, by 
considering defects in an n-fold degenerate ground state, 
analogously to the generalization of the polymer system 
0, EH- To illustrate, we consider the fractional fermion 
number 1/3. We consider the ground state with a 1/3 
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FIG. 7: The nonuniform spatial profile for the hopping am- 
plitudes, (a) Here '±' denote the couplings a ± /i in Eq. (|6]) 
where one of the signs of the hopping amplitudes has been 
changed. Introducing the sign switch generates two defects 
that are marked by the vertical lines, (b) Here '±' denote 
the coupfings a ± /i in Eq. (|42} , with A^^' given by Eq. (|43j , 
where one of the signs in the middle has been switched. In 
this case the sign change introduces three defects (marked 
by vertical lines). The first defect represents the boundary 
between the p = and p = 2 three-fold degenerate configura- 
tions of Eq. (|43|l : the second between p = 2 and p = 1; and 
the third between p = 1 and p = 0. Part (c) represents the 
same configuration as shown in (b), but with the three defects 
moved apart. 



atomic filling factor. We assume that the hopping am- 
plitude between the lattice sites reads 



(42) 



Here the phase factors A'^^ exhibit a periodicity of three 
lattice sites and define a three-fold degenerate ground 
state (in the limit of an infinite lattice) for the values 
p = 0, 1, 2. We choose 



A 



(p) 



1 



47r(/c -I- p) 



3 3 

Then for p = we have 

x + x + x — x + x + x — x + x + x — X 



(43) 



(44) 



where we again use an x to denote a lattice site and ± 
the couplings a ± ^. Similarly, for p ~ I we obtain 



x + x — x + x + x — x + x + x — x + x- 



(45) 



In the dimerized lattice system with the coupling ([5]), 
one may easily observe that changing the sign of one of 
the hopping amplitudes corresponds to adding two identi- 
cal defects in the systems [Fig. [7l[a)] and by changing the 
number of atomic states in the vacuum by one (for the 
ground state of the average of 1/2 atoms per site). The 
half fermion results from sharing the extra state evenly 
between the defects. Similarly, one may show that chang- 
ing one of the signs of the hopping amplitudes in Eq. ((^ 
corresponds to adding one atomic state in the 1/3 filled 
vacuum. It is relatively easy to see that such a sign switch 
can be viewed as three identical defects that join together 
the three degenerate vacuum configurations of Eq. (|43p . 
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representing the boundaries between the p = and p = 2, 
p = 2 and p = 1, as well as the p =1 and p = states; see 
Fig.jT^b-c). If the defects are moved wide apart compared 
to the size of the kink bound states, the extra atomic 
state has to be split equally between the defects, result- 
ing in the 1/3 fcrmion number. 



IV. LOW ENERGY CONTINUUM LIMIT AND 
FRACTIONALIZATION IN RELATIVISTIC 
FIELD THEORY 

In the continuum limit the Hamiltonian ([1]) [with 
r^(r) defined in Eqs. p7|) and ([TS]) ] can be transformed 
[45l.l46l. into a relativistic Dirac Hamiltonian exhibit- 
ing fractional fcrmion number in quantum field theory 
0. The transformation is mathematically the same as 
in the ID polymer case 0, [H, m, H 1131 , but the physics 
is rather different: in our atomic system the kink (p{x) 
appears in the em-induced coupling between the internal 
atomic states, and not as a physical domain wall kink. 

A particular advantage of the continuum quantum field 
theory is that it is amenable to a simple description. The 
continuum limit of the low-energy expansion corresponds 
to the linearization of the band structure about the Fermi 
surface by writing e(p) ~ e{pF) + {p~PF)^' {pf), with the 
Fermi momentum pf = 7r/2(i, and keeping only terms to 
leading order in small d. It becomes accurate in the dilute 
gas limit, when the atomic correlation length ^ ~ \ad/ 
from Eq. (|15p . is much larger than the lattice spacing. 
In the continuum limit we write the fermionic annihila- 
tion operators for the even and odd sites as continuous 
functions of the lattice spacing d: 

C2fc+i = (-l)(2''^+i'/2VM7/[(2fc -I- l)d], (46) 
C2k^{-lfV2dv{2kd). (47) 

With these identifications, the continuum limit proceeds 
exactly as in the polymer case H. To leading order in 
small d we obtain 

u'{2nd) ~ [u{2nd + d) - u{2nd - d)] /2d, (48) 

v'i2nd + d)~ [v{2nd + 2d) - v{2nd)] /2d, (49) 

u{2nd) ~ [u{2nd + d) + u{2nd - d)] /2, (50) 

v{2nd + d)~ [v{2nd + 2d) + v{2nd)] /2 . (51) 

Here u'{nd) denotes a discrete spatial derivative of w. By 
setting Ck = {—l)''S/2 and with the hopping determined 
by Eq. ([TT]) . the Hamiltonian ([1]) reads 

H /h = -2id'^a^[u'' {nd)v' (nd) + v\nd)u' (nd)] 

n 

Sd 

+ — [u^ {nd)u(nd) — {nd)v{nd)] 

n 

+ 2i^d (p{nd) [u^ (nd)v(nd) — v\nd)u{nd)] . 



In the continuum limit we replace nd —f x and dX^n ^ 
J dx. Combining the two states into a two-component 
spinor 



(53) 



(52) 



and making the transformation: ^ — > exp (i7r(T'^/4)^', 
we may express Eq. ([5^ as the relativistic (H-l)- 
dimensional Dirac Hamiltonian 0,0] for the spinor ^{x) 
coupled to a bosonic condensate 'p(x), 

H= I dx[cfi^^a'^^ + hgip-^'^a^-^+mcH^aH]. (54) 
J idx 

Here we have identified c = 2da, g = 2/i, and m = 
hS /{Sd'^a'^). The cr* denote the Pauli spin matrices, g 
the coupling coefficient, m the fermionic mass, and c the 
speed of light. In the atomic system the spinor compo- 
nents refer to the two internal atomic levels. The corre- 
sponding eigenvalue system reads 



ch— h hgip ) u{x) — mc^v{x) =Ev{x), 



— cfi.— + hgip ) v{x) + mc^u{x) =Eu{x) . (55) 



Note that in the relativistic Dirac equation ([55)1 for the 
fermionic atoms both the speed of light and the coupling 
coefficient may be controlled by the different em field am- 
plitudes a and fi. The mass of the Dirac fermions vanishes 
at the exact resonance coupling between the two atomic 
levels and is linearly proportional to the detuning that 
can be accurately changed in the experiments. More- 
over, in the relativistic theory, as in our atomic scheme, 
the bosonic field f{x) can be taken to be a static classical 
background field. Thus, the atomic Hamiltonian ([1]) with 
a defect in the pattern of alternating hopping amplitudes 
[Eq. (|18p ] can also be written in the continuum limit in 
the form jM]) and ((55)) . 

Fractionalization in the relativistic theory is described 
as follows [2, 1^ . Suppose the bosonic field (p in Eq. ([51]) 
has a doubly-degenerate ground state, for example aris- 
ing from a double-well potential V{(p) oc {ip-^ — 7^)^. 
There are two minima, f{x) = ±7, and the reflection 
symmetry (f —ip is spontaneously broken. Conse- 
quently there are topological solitons (p{x) (also known 
as "kinks" ) interpolating between these two degener- 
ate minima. For weak coupling, we make a Born- 
Oppenhcimer approximation, so that the fermion number 
must now be defined through the second quantization 
of the fermion field in the presence of this background 
soliton. Being a Dirac particle, the fermion has both 
positive and negative energy single-particle eigenstates, 
and therefore an associated Dirac negative energy sea 
in which the negative energy states are occupied. This 
Dirac sea is the ultimate physical origin of the fractional 
fermion number [4|. 
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The second-quantized number operator can be defined 
in two equivalent ways. The formal relativistic quantum 
field theory definition is [i^ 



(56) 



with the commutator ensuring that the fermion particle 
number vanishes in the free vacuum. Equivalently, one 
defines the physical particle number in the soliton sector 
as being measured relative to the free vacuum sector, so 
that the fermion number density is 



p{x) 



dE{\^E{x)\'-\M^)\') 



(57) 



The same conclusion is reached Q using the subtracted 
definition (j57p of the number density operator, if we sim- 
ply combine the expression (|57p with the completeness 
of the two sets of eigenstates, {^'b} and {i^e}, h^ the 
soliton and vacuum sector, respectively: 



/O r+oo 
d£;|*B(x)|2 + |*o(a;)|'+ / dE\^E{x)\ 
-oo J0+ 



-t-oo 



dE\^Eix)\' 



(61) 



Thus, in the fermion-soliton vacuum defined with only 
negative energy states filled (i.e., with the zero mode 
empty), the fermion number density (|57|) is 



where ^e {"^e) are the fermion single-particle energy 
eigenstates in the soliton (free) vacuum sector. Once 
again, by construction the fermion particle number van- 
ishes in the free vacuum. 

It is easiest to exhibit the fermion number fractional- 
ization when m = in the Dirac Hamiltonian (|54p . In 
this case the Hamiltonian has a conjugation symmetry 
that pairs positive and negative energy states. For every 
eigenvalue e of Eq. ([55]) there exists an eigenvalue — e, 
and the corresponding eigenfunctions are paired accord- 
ing to -c = cr^'^l- However, the nontrivial feature of 
this fermion-soliton system is that there is also a zero- 
energy bound state ^0(2;), and furthermore, this state is 
localized at the soliton jump: 



'J'o(.t) 



A exp [- 



(58) 



This state is self-conjugate, cr'^^o = ^0, and results in a 
doubly-degenerate soliton sector vacuum 0]. 

Thus, the mode expansion of the second quantized 
fermion field operator has the form 



*(2;,t) = a*o(a;) 



(59) 



dk 



where hk and Ck denote annihilation operators for contin- 
uum fermion and antifermion modes (respectively), while 
the operator a is the annihilation operator for the self- 
conjugate zero-energy state. Then the fermion number 
operator ([551) ™ the presence of the soliton reads Q: 



N^\[a\a] + \ j dk {[blM + [ck.cl]) 



+ / dk{blbk~clck), 



(60) 



Note that the operators a and couple the two degen- 
erate fermion-soliton ground states in which the zero en- 
ergy single-particle state is either occupied or not. Thus, 
it follows from (|SD)) that the ground-state fermion-soliton 
states possess fractional fermion numbers ±1/2. 



p{x) 



-2l*o(:c)P 



(62) 



which yields fermion number —1/2. Similarly, if the zero 
energy state is filled, the fermion number is -1-1/2, con- 
sistent with the result above from ([60]) . 

In fact, it is not just the expectation value of the 
fermion number that is half-integer; with a suitable 
physical definition of the number density operator as a 
smeared operator defined over a physical sampling re- 
gion, it can also be shown that these half-integer values 
are eigenvalues of the number operator [i^ . [4^. The 
fermion-soliton system yields eigenstates of this physical 
number operator. Furthermore, the fractional part of the 
fermion number has a topological character: it is insen- 
sitive to local deformations of the bosonic soliton field, 
depending only on its asymptotic behavior Q . 

For m 7^ 0, the situation becomes even more interest- 
ing [l3| . In this case, the conjugation symmetry of the 
Dirac Hamiltonian (|54p is broken and the positive and 
negative energy states are no longer paired in a simple 
way. In the limit of a slowly varying soliton (on the scale 
of the fermion Compton wavelength), the soliton acts as 
an inhomogeneous electric field that polarizes the Dirac 
sea vacuum, thereby building up local fermion number 
at the location of the kink. A straightforward one-loop 
computation [l3, 49, 50j] yields a fermion number taking 
arbitrary fractional values: 



1 / ^97 

N ^ arctan 



(63) 



In our atomic system studied in Section IIB, we studied 
the conjugation symmetric case where = in Eq. (j24p . 
which corresponds to the conjugation symmetric m ~ 
case in the continuum limit. If instead we set a nonvan- 
ishing effective detuning S for the hopping field between 
the two atomic levels, so that = {—1)^6/2 with (5^0, 
the symmetry between the valence and the conduction 
band states is broken, and they are no longer coupled 
together in a simple way. In the continuum limit this 
corresponds to the m ^ case in which the conjuga- 
tion symmetry is broken. Then the fractional part of the 
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fcrmion number may exhibit any value, in the continuum 
limit according to 



1 /4/i 

A* = arctan — 

TT \ 



(64) 



The ratio 4/i/(5 between the coupling strength and the ef- 
fective detuning therefore determines the fractional part 
of the particle number. In experiments this could be 
engineered accurately, allowing a potentially controlled 
way of preparing the fractional part of the eigenvalues 
by tuning the mass of the Dirac fermions. 

An interesting further variation would be to study the 
temperature dependence of the fermion number. The 
above discussion is all for zero temperature, but at 
nonzero temperature the fermion number may develop 
nontopological contributio ns [14| and the fluctuations 
may become nonvanishing [15|, depending on the form 
of the kink. Nonetheless, the fractional fermion is not 
a singular feature of the theory at T = that vanishes 
without a trace at any nonzero temperature. As a con- 
crete example, for a kink the finite temperature induced 
fermion number is [ij, [IE, HH, HI, HI] 



sin Q 



^ ,1^ ((2?^ + 1)' + y') ^(2?^ + 1)' cos2 + y2 



where Q = arctan(4/^/(5), and y = mc^/(7rT). This ex- 
pression smoothly reduces to the zero temperature result 
dSll) as T ^ 0. The fluctuations (AiV)^ = {N'^) -(N)'^ 
also has a characteristic temperature dependence |15| . 

The crucial part of our proposal for fractional particle 
number is the em field '^(x) in Eq. (|18p . This is very dif- 
ferent from the fermion particle number fractionalization 
in polymers, as our fermionic fields are not coupled to a 
bosonic matter field with a domain-wall soliton. Instead, 
the coherent em field with a topologically appropriate 
phase profile is coupled to the FD atoms via internal 
transitions. This results in the quantization of the FD 
atomic gas with nontrivial topological quantum numbers 
corresponding to the soliton sector of the relativistic 1-1-1 
quantum field theory models of fractionalization. On the 
other hand, a spatially constant phase profile ^p{x) rep- 
resents the FD vacuum sector exhibiting integer particle 
numbers. 



V. OPTICAL DETECTION 

The dimerized optical lattice resulting from the alter- 
nating pattern of the hopping matrix elements causes the 
single-particle density of states to acquire an energy gap, 
which in the limit iV^ — > oo equals 4fi|/i|; see Eq. ([H]). 
The zero state is located at the center of the gap. The 
resulting excitations at half the gap energy could be de- 
tected by resonance spectroscopy [5J]. This provides in- 
direct evidence of fractionalization, as in the polymer 
systems [l^. Because in our scheme [l| the gap is pro- 
portional to the amplitude of the em field inducing the 



hopping, the size of the energy gap can be controlled ex- 
perimentally. 

Such mid-gap spectroscopy, however, would not pro- 
vide information about fermion numbers. In this section 
we consider direct optical measurements of the expecta- 
tion value and of the fluctuations of the fermion number 
in order to ascertain if they are compatible with fraction- 
alization. 

We assume that far off-resonant light excites the 
atoms, and consider the ID optical lattice to be opti- 
cally thin. We take the light scatterer at each lattice 
site to be much smaller than the wavelength of the de- 
tection light. In the case of off- resonant excitation the 
amplitude of the light scattered from an essentially point 
source is proportional to the number of atoms [2^. Also, 
for off-resonant excitation the scattering is coherent; the 
scattered light has a fixed phase relation to the incom- 
ing light, and is fully capable of interference. Thus, in a 
given point of observation the positive frequency parts of 
the field operator for the light scattered from the lattice 
sites simply sum up to 



(65) 



C is a constant containing the overall intensity scale of 
the detection light. The factors ak include aspects such 
as intensity and phase profiles of the detection light, ef- 
fects of the spin state at each site k on light-atom cou- 
pling, and changes of the amplitude and phase of the 
light as it propagates from the lattice sites to the point 
of observation. A more detailed description how Eq. (pS]) 
is obtained can be foimd in Ref. 1251. 



A. Measuring atom numbers 

In forward scattering and variations thereof such as 
phase contrast imaging, the scattered and the incoming 
light interfere. The ultimate measurement of the inten- 
sity in effect records the expectation value of the scat- 
tered electric field E = (E^). The observable at the 
detector is 



E 



k 



ak{c\ck) 



kp 



akUf^pUp. 



(66) 



This is a linear combination of the expectation values 
of the numbers of fermions at each lattice site with the 
coefficients ak, exactly as introduced in Eq. ([M)) . 

We now construct a numerical example about forward 
scattering as a means to detect the fractional fermion. 
We make use of the fact that the fermion species at the 
alternating lattice sites are likely to be different. We 
assume that the detection light is far blue-detuned in one 
species and far red-detuned in the other, and that the two 
dipole matrix elements are comparable. One may then 
find a laser tuning such that the intensity of the scattered 
light is the same for both species. However, the lights 
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scattered by the two species are out of phase by tt, and 
out of phase with the incident hght by ±7r/2. With the 
usual tricks of phase contrast imaging, the relative phase 
of incident and scattered light is then adjusted so that 
in interference the light from one species directly adds to 
the incident light, and the light from the other species 
subtracts. We incorporate the alternating sign into the 
definition of our observable, and write 

i? = ^a,(-l)^4cfc). (67) 

k 

We take it that the incoming light is a plane wave with 
a constant phase and amplitude. It will therefore just 
contribute a common constant into the envelope factors 
Q!fc, hence basically stands for the amplitude of the 
light as transmitted through the imaging system from 
the source point, lattice site r^, to the observation point, 
r. We adopt a rudimentary physical-optics model accord- 
ing to which an imaging lense first takes a Fourier trans- 
form of the light field at the object plane, passes only the 
Fourier components that make it through the aperture 
of the lense, and then takes another Fourier transform to 
form the image. 

Let us say that the geometry has been arranged in 
such a way that all Fourier components of light in the 
plane of the aperture up to the absolute value K are 
passed, the rest are blocked. Depending on where the 
object and image planes are, there might be a scaling of 
the image with respect to the source, but we ignore both 
this and the usual inversion of the image with respect to 
the object. The transfer function of the lense from point 
r/. to the point r, normalized to the peak value of unity, 
is then 2Ji{K\r — rk\)/{K\r — rfc|), and so we have the 
transmitted field 

We present an example in Fig. [H] We choose the 
parameters fi ~ —0.1a, and the numbers of sites and 
fermions Ng = 129 and Nf = 65. Given the imaging 
light with the wavelength A, we use the cutoff wave num- 
ber K = 2n/{V5X). This would be appropriate, for in- 
stance, if the imaging system had the numerical aper- 
ture F = 1 and the lattice resided near the focal plane 
of the imaging lense; or if the numerical aperture were 
the rather extreme F = 0.5, and both the optical lattice 
and its image were removed by two focal lengths from 
the lense. We assume that the wavelength of the lat- 
tice light and of the imaging light are the same and that 
the spacing between the lattice sites is d = A/4. Under 
these circumstances it would difficult to resolve individ- 
ual lattice sites. In practice a ID optical lattice can be 
significantly stretched by changing the intersection an- 
gle between the laser beams, and so the resolution of the 
optical detection is not necessarily limited by the lattice 
site spacing even if the wavelengths of the lattice light 
and the light used for optical detection were comparable. 
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FIG. 8: Optical image (dashed line) of the half-fermion den- 
sity profile on top of a uniform background density with the 
average number of 0.5 atoms per lattice site produced by a 
specific imaging system, as discussed in the text. The density 
profile (solid line) that is being imaged is also shown. The 
size of the soliton is set by the choices Na = 129, Nf = 65, 
and n = —0.1a, corresponding to the correlation length of 
5 ~ 10 lattice sites. 

We set c? = A/4 simply to limit the number of parame- 
ters to consider. We plot the optically imaged fcrmion 
lattice along the line of the atoms (dashed line), and the 
number of fermionic atoms in excess of the average occu- 
pation number i for the even-numbered sites that carry 
the lump with (N) = ^ (solid line) as obtained from 
Eq. (|34p . The curves are normalized so that the maxima 
overlap. 

Although the quantum operator describing the 
phase contrast imaging is not directly proportional to 
the fermion number operator TV defined by Eq. (|34p . 
its expectation value produced by the imaging system, 
Eq. nevertheless accurately depicts a resolution 

rounded version of the half-fermion lump. In fact, phase 
contrast imaging has been used for nondestructive mon- 
itoring of a Bose-Einstein condensate (ssj . and the ab- 
sorption of a single trapped ion has been detected exper- 
imentally long ago [5^ . While a lot of assumptions went 
into our specific example, a light scattering experiment 
along these lines should be feasible with the technology 
available today. 

B. Measuring atom number fluctuations 

1. Intensity of Bragg peaks 

In the absence of interference with the incoming light, 
the straightforward observable is the intensity of the light 
scattered from the atoms. We have 

I^{E-E+). (69) 

With illumination of the optical lattice by a focused 
laser beam it is possible to arrange the field strength of 
the detection light vary from site to site approximately 
as a Gaussian. Let us further assume, contrary to the 



previous section, that the atoms at the lattice sites are 
identical as far as light scattering is considered. If now 
the detection is carried out in the direction of construc- 
tive interference (Bragg scattering) so that the light am- 
plitudes scattered from the lattice sites have the same 
phase, and far enough from the lattice so that the prop- 
agation distance of light from each lattice site to the de- 
tector is approximately the same, the field strength is of 
the form 

E+ =Y^e-^^/-^\lck. (70) 

k 

A measurement of the intensity of the light therefore 
measures the square of the atom number operator, / = 
(TVtAT) = |(iV)|2 + {/^Nf, as in Eq. with the Gaus- 
sian weights ajf. 

The key point in the detection of the correlated atom 
number fluctuations responsible for the half-fcrmion is 
to measure the coherent fluctuations (|37p . i.e., to rely on 
interference of light scattered from different lattice sites. 
If the light scattered from individual sites can be resolved 
or if a too broad angular average in the detection or other 
such cause wipes out the interferences [a^aj 5kial.OLk], 
the scattered light probes the "incoherent" fluctuations 
{/S.N)i [Eq. ([^T]) ]. as if the fermion number fluctuations 
in each site were independent. 

The problem in detecting the fractional fermion num- 
ber is that in the measured signal of the light intensity the 
fluctuations AA^ appear on a nonzero background. We 
demonstrate by plotting in Fig. [5] separately the contri- 
bution |(A^)p, as if the fermion numbers were precisely 
fixed, and the fluctuation term (AiV)^. We also show 
the fluctuations (AA^)^ that would result if the fermion 
number fluctuations at adjacent sites were uncorrelated. 
These are given as functions of the width of the focus 
w of the laser beam. Here Ng = 129, Nf = 65, and we 
choose n = —0.9 a to make a sharply localized zero state. 

In our discussion we employ the Gaussian width w = 4, 
which would ordinarily represent very tight focusing. 
Here the contribution from atom number fluctuations 
to light intensity is two orders of magnitude below the 
coherent intensity, whereas the fluctuations from uncor- 
related fermion numbers would make a contribution an 
order magnitude smaller than the coherent intensity. As 
our detection light was assumed to be far off resonance, 
photon number fluctuations are Poissonian. Under other- 
wise ideal conditions, including absolute knowledge and 
control of all experimental parameters pertaining to in- 
tensity of the scattered light, the detection of about a 
hundred photons could reveal the difference between cor- 
related and uncorrelated fermion numbers, whereas a 
quantitative study of the actual correlated fermion num- 
ber requires the detection of about 10,000 photons. Un- 
fortunately, a large number of scattered photons means 
a large number of recoil kicks on the fermions. Currently 
available optical lattices likely cannot absorb the assault 
of hundreds of photon recoils without developing some 
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FIG. 9: The intensity of light scattered from the optical lattice 
if the fermion numbers did not fluctuate, |(7V)|'^, and the ad- 
ditional intensity due to fermion number fluctuations (AA*')^, 
as a function of the size of the focus w of the driving light 
given in the lattice units. We also display the added inten- 
sity (AA)f that would result if fermion number fluctuations 
were uncorrelated between adjacent lattice sites. The soliton 
parameters are Ns = 129, Nf — 65, fi = —0.9 a. 

form of a dynamics that complicates the phenomena we 
are analyzing. 



2. Variations of intensity measurements 

While the scattered light in principle conveys informa- 
tion about the fluctuations of atom number, the method 
to extract the information we have discussed so far is 
quite challenging. 

A way to discard the uniform background of ^ 
fermions per lattice site, as in going from operator N 
to A^ in Eq. would be most useful. One might, for 

instance, think of the alternating-sign Gaussian envelope 
a^, as defined in Eq. (j40[) and (in principle) realized as 
discussed in Sec. IV A| and still look in the direction of 
ordinary Bragg scattering. By virtue of the first form of 
Eq. the coherent signal |(A)P from the background 
of i indeed cancels in the limit of a broad envelope. How- 
ever, the fiuctuation part {AN)^ as appropriate for is 
not small, and does not directly bear witness to the small 
fluctuations of the atom number of the half-fermion. On 
the contrary, as discussed in Sec. Ill B 21 in this case the 
fluctuations of the atom number under the Gaussian en- 
velope are larger that one would expect if the atom num- 
ber fluctuations at adjacent sites were uncorrelated. 

A scheme that fares better is based on the observation 
that off-resonant light scattering is coherent. If the atom 
numbers were fixed, the field scattered by the lattice as 
a whole were also completely coherent with the incoming 
light. One could then tap part of the incident light and 
impose a suitable attenuation and phase on it so that on 
the detector the incident light exactly cancels the scat- 
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tered light. Given the fiuetuating atom numbers, it is 
correspondingly possible in principle to arrange things so 
that the detected intensity originates entirely from atom 
number fluctuations. 

In technical terms, assume that the incident and scat- 
tered fields are superimposed on the detector so that the 
electric field operator is 

E+ = cY,ak{Ack-{clck)) (71) 

k 

instead of Eq. (pS)) . This eliminates the coherent part of 
the intensity cx |(7V)p, leaving only the fluctuation part 
(AA^)2 in Eq. The situation of Eq. dZH) may be 

realized operationally by minimizing the detected inten- 
sity with adjustments of the attenuation and the phase 
of the canceling light. The remaining intensity is then 
directly proportional to the the square of atom number 
fluctuations under the envelope of the detection light. A 
reduction in the minimum detected intensity with an in- 
creasing size of the focus of the driving light would be a 
signature of anomalous atom number fluctuations char- 
acteristic of the fractional fermion. Another approach to 
eliminate the coherent part of the light intensity is to ex- 
ploit the diffraction pattern of the scattered light from 
the regular array of the lattice sites and detect the pho- 
ton number fluctuations at the location of the destructive 
interference of the coherent cx contribution of the 

scattered light. 

Although our goal here is not a specific experimental 
design, another variation with potential to overcome the 
accumulation of atom recoil merits a mention. So far 
we have dealt with what in essence is spontaneous Bragg 
scattering. Recently, induced Bragg scattering has been 
introduced as a method to study the condensates in de- 
tail [13, nil . In optical lattices the strong spatial confine- 
ment r riay complicate the Bragg spectroscopy measure- 
ments [1^, [gO] , but the advantage of the light-stimulated 
Bragg scattering is that conceivably the light pulses could 
be made so short that the harmful effects of photon recoil 
do not have time to build up during the measurement. 



Hamiltonian ^ . The length scale over which the energy 
variation due to the trap becomes comparable to the en- 
ergy gap is then = (8?i|/x|/Mw^)^/^, where lo denotes 
the trap frequency and M the mass of the atom. For 
the system to have a locally homogeneous gap value, this 
length scale should be much larger than the correlation 
length ^. Moreover, in order to have a better experimen- 
tal control over the atom numbers and to maintain the 
one-half filling throughout the lattice, we require that 
the energy gap is nonvanishing everywhere, or tq > Nhd. 
The nonvanishing gap parameter and the locally homoge- 
neous limit could therefore be reached with a sufficiently 
weak trap and by experimentally controlling the size of 
the gap. By using sufficiently large values of it might 
also be possible to study the effective zero temperature 
limit ^ ksT/h and to have a better control over the 
light scattering measurements. Finally, the sharp edges 
of the optical lattice system with hard-wall boundaries 
could be prepared, e.g., by shining blue-detuned laser 
beams at the lattice ends. 

We could possibly also construct complex higher di- 
mensional models in 2D or 3D optical lattices with atoms 
coupled to nontrivial em fields. In relativistic (2+l)D 
quantum field theory, a fermionic field coupled to a 
bosonic field exhibiting a vortex profile results in fermion 
particle number fractionalization [3, [6ll, [l^, [H, [G^I • It 
may also be possible to find analogies to relativistic 
(3-|-l)D quantum field theoretical models 0, e.g., by us- 
ing em field superpositions to prepare topologically non- 
trivial atomic configurations [s^, |43, H^, E^, [63, [11] . The 
significance of localized zero energy modes in fermionic 
systems has also recently been discussed in the context of 
quantum teleportation [69|, usin g re sults from the quan- 
tum field theoretical analysis in |70 |. 
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